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1. Introduction 

Since overlap fermions have an exact chiral symmetry on the lattice, they are an appropriate 
tool to investigate the chiral phase transition.^ Simulations involving the effects of dynamical 
overlap sea quarks are just becoming possible now, but as they are still numerically extremely 
demanding on realistic lattices, we use a hybrid approach by implementing the overlap action for 
the valence quarks on dynamical configurations generated with Nf = 2 flavours of ff{a) improved 
Wilson sea quarks. 

As the lattice spacing, a, for dynamical fermions is connected to both the coupling constant j8 
and Ksea, increasing the temperature T = l/{N4a) across the transition is affected by an increase 
in Ksea at fixed values of the coupling, j8, and time-extent, N4, of the lattice. In order to probe the 
topological properties of dynamical gauge fields in the vicinity of the phase transition, we con- 
centrate on the values Ksea = 0.1343, 0.1348 and 0.1360 corresponding to T/Tc = 0.98, 1.06 and 
1.27, respectively, at fixed j8 = 5.2 on 16^ x 8 lattices generated by the DIK-collaboration [§]. The 
lattice spacing is approximately <2 = 0.12 fm. We use (^(200) configurations for each value of 
Ksea- The transition point Kc can be identified as the point where the Polyakov loop susceptibility 
X = Ns{{L^) — {L)^), with the Polyakov loop L(s) = ^TrH^^i C/4(?,^4) , assumes its maximum. 
Recent work of the DIK-collaboration on 16^ x 8 lattices at j8 =5.2 gave = 212(2) MeV, corre- 
sponding to Kc = 0.13444(6) [|]. 

2. Density and locality of the eigenmodes of the overlap operator 

The spontaneous breaking of chiral symmetry by the dynamical creation of a nonvanishing 
chiral condensate, (^^), is related to the spectral density p(A) of the Dirac operator near zero 
by the Banks-Casher relation (^^) = — ^p(O). Using the Amoldi-algorithm, we compute the 50 
lowest eigenvalues on every configuration.^ To smooth the configurations we perform 1 step of 
APE-smearing, using the smearing coefficient a = 0.45. In Fig. |l| we plot the spectral density. 
One can clearly see a nonvanishing chiral condensate in the chirally broken phase at Ksea = 0.1343 
below the phase transition and a large gap in the spectrum at Ksea — 0. 1360 in the symmetry restored 
phase. However, at Ksea = 0.1348, a value which was determined by the Polyakov-loop method to 
be above the transition, we find a nonvanishing tail of small eigenvalues extending down to zero. 

A useful measure to quantify the localisation of eigenmodes ^ is the inverse participation 
ratio / = VJ^^p(x)^, with the scalar density p{x) = wl{^)¥x(^) for normalised eigenfunctions 
Y,xP{^) = 1- While / = y if the density has support only on one lattice point, / decreases as 
the density becomes more delocalised, reaching 7=1 when the density is maximally spread on 
all lattice sites. Fig. ^ shows the IPR's for all configurations as a function of the modulus of the 
eigenvalue below and above the phase transition. One can see that in both cases the zero modes and 
the first non-zero eigenmodes are highly localised, while the mean IPR for the highest computed 
eigenvalues is approximately 1.3. 

^ Our results on the vacuum structure for T = quenched configurations are presented in the talk by Y. Koma [|l|] at 
this conference. 

^For our implementation of the overlap operator see Ref. [Q. Also here we use p = lA and stereographically 
project the eigenvalues. 
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Figure 1: The spectral density of the overlap operator in the vicinity of the phase transition (50 lowest 
eigenvalues included). 
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Figure 2: The IPR of overlap eigenmodes (50 eigenmodes from all configurations) below and above the 
transition. 

3. The local structure of topological charge fluctuations 

The topological charge density for any 75-Hermitean Dirac operator satisfying the Ginsparg- 
Wilson relation is defined as ^\ 

q{x) = ^Tr ysD{x,x), Q = ^^(x). (3.1) 
To compute the topological charge density, we use two different approaches [0, |]. In the first 



approach, we directly calculate the trace of the overlap operator according to equation Eq. ([3.1|). 
This is a computationally very demanding task and is therefore performed on only 5 configurations 
for each Ksea- The "full" density computed in this way includes charge fluctuations at all scales. The 
second technique involves the computation of the topological charge density based on <^(50) low 
lying modes of the overlap Dirac operator. Using the spectral representation of the Dirac operator, 
the eigenmode expansion of the topological charge density reads 

= -1^(1 - \) c^x), c\x) = yifl{x)ysWx{x) • (3.2) 

Truncating the expansion acts as an ultraviolet filter by removing the short-distance fluctu- 
ations from q{x). Note that the topological charge Q = Y.xQx{^) is not affected by the level of 
truncation. 
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Figure 3: Isosurfaces of topological charge density with \q{x)\/qmax = 0.1, 0.2, 0.3, 0.4, with red (green) 
surfaces indicating positive (negative) charges in one timeslice of a configuration with Ksea = 0.1343. The 
upper (lower) pictures are based on the full (eigenmode truncated) density. 



Reflection positivity, i.e. positivity of the metric in Hilbert space, demands that the topological 
charge correlator is negative 

C<iir) = ^Y.iq{x)q{y))<0, r>0, (3.3) 



X 



where r — \x — y\ is the Euclidean distance. Since overlap fermions are not ultralocal, the fermion 
action is not strictly reflection positive. Therefore, analysing dynamical or quenched lattices with 
valence overlap fermions, one may expect a positive core and a negative tail of Cq{r) [|l|, |l^]. This 
will be washed out when a finite UV filter is applied. 

Let us first consider the local structure of the topological charge fluctuations, studying sign co- 
herent clusters formed by connected neighbours Xi with \q{xi) \ > qcut = 0.1 . . .0.4^^^;^, where qcut 
is a cutoff value to be varied at will. To give a picture of the spatial distribution of the topological 
charge density, isosurface plots with q{x) = -^qcut are shown in Fig. ||.^ 

It is obvious that the UV cutoff in the topological density results in a completely different 
structure actually seen by the lowest modes. This is also reflected in Fig. ^ which shows the 



topological charge correlator Eq. ( |3.3| ) in terms of the truncated density. We take 50 modes into 
account, which seems not to be sufficient to generate the negative tail of the correlator [jlj]. Since 
hypercubic rotation symmetry could be broken we plot the correlator as a function of the Euclidean 
distance r and the temporal distance t. The UV truncated correlator represents mainly size, shape 
and temperature dependence of the clusters seen in the lower row of Fig. ||. We see a slightly more 
negative correlation developing in the high-temperature phase. 

Further investigating the full density, we show in Fig. ^ (a) and (b) that the number of clusters 
reaches a maximum at a moderate cutoff value, while for a cutoff as low as qcut =0.1 qmax there 



Movies are available at [http : / /www. cip .physik . uni-muenchen . de/ "weinberg/topdens / ] 
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Figure 4: The topological charge correlator using 50 eigenmodes on ^(200) configurations. 
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Figure 5: The total number of clusters of full density (a), the fractional size of the largest cluster (b) and the 
distance between the 2 leading clusters (c) depending on the cutoff value. 



are 2 totally dominating oppositely charged clusters with a fractional volume of around 50% each.^ 
The fractional volume and the total charge of the clusters ordered by their size are plotted in Fig. ^. 
Studying the point-to-point correlation of the largest clusters based on the full charge density, we 
find that the clusters start to percolate at qcut = 0.2 (0.25) qmax for Ksea = 0. 1343 (0. 1360). Fig. ^(c) 
shows the maximum of the minimal distance between each point of the largest cluster and every 
point of the second largest one. The small value for low cutoffs in q implies that the 2 dominating 
clusters are tangled and intertwined in a complex way. 

A preliminary dimensional analysis indicates that we are rediscovering for the full density the 



picture of 2 dominating approximately 3D sign coherent clusters QllD at low cutoffs also at finite 
r, whereas we see ID highly charged small clusters at high thresholds. This picture means that 
the QCD vacuum model with 4D coherent (anti)instantons with a typical instanton radius of 0.3 — 
0.4 fm is strongly modified by quantum fluctuations. The eigenmode-truncated density is more 
compatible with the conventional multiple-lump picture. Apart from the topological susceptibilty 
which is strongly suppressed in the chiral- symmetry restored phase, only quantitative structural 
differences between the phases are seen so far. 



^The maximal number of clusters of the truncated density is an order of magnitude smaller. 
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Figure 6: Fractional volume (uppper plot) and total charge (lower plot) of the clusters of full charge sorted 
by their size for qcut = 0.10 qmax (left) and 0.40 qmax (right) for one configuration with Ksea = 0.1343. 
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